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1 Introduction 



Systems of partial differential equations of the form 

ri, y = e^^ rj , i 1.2 V. (1.1) 

called generalized two-dimensional Toda lattices have very important applications in Liouville and 
conformal field theories, they are studied in details (see [1]-[T2] and the references therein). Here 
the matrix A = {dij} is the Cartan matrix of an arbitrary finite or affine Lie algebra. It is known 



that in the former case system (1.1) is Darboux integrable while in the latter case - S-integrable. 



The widely known Drinfel'd-Sokolov formalism allows one to construct the Lax representation for 



the system (1.1) in terms of the Lie algebra canonically associated with the corresponding Cartan 
matrix A. 

The problem of finding discrete versions of the system is intensively studied (see, for instance, 



13]- [26]). Recently, in [27] integrable differential-difference analog of system (1.1) was suggested 



v [ x - 4 = e E ^ 1 a ^ vi +^j=i+^A+hi,i^+v\) ^ i = i j2 , N. (1.2) 

Here the functions u* = v^(n,x), j = 1,...,N are the searched field variables. The subindex 
denotes a shift of the discrete variable n or the derivative with respect to x: r 3 k = r^(n + k, x) and 
r J x = ^r^{n,x). A new example of the semi-discrete system corresponding to Cn was recently 
studied in 



In the present article we study the problem of further discretization of system (1.1), i.e. the 



problem of finding a rule allowing to assign to any Cartan matrix a system of integrable difference- 



difference equations approximating in the continuum limit system (1.2) and therefore system (1.1). 
As it was pointed out in [T3] the problem of discretization is important from physical viewpoint, 
they might have applications in discrete field theory and in quantum physics (see also |15]-|17j). 
They can also be regarded as difference schemes in numerical computations (see [31]). 

Following questions were addressed in |14j : 
1) whether there exists an integrable discrete version for any two dimensional Toda field equation 

0); 



2) which kind of algebraic structure (like Lie algebra or Lie group) is naturally related to discrete 
versions. 
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The subject has intensively been studied during the last 10-15 years. Various discrete versions 
of Toda field equations were investigated in literature (see [16], [17], [18], [19], [21] and references 
therein). As an alternative answer to the first question we suggest a fully discrete version of system 
( 1.1[ ) in the following form: 

e -<l+ u i,0+<l- u 0,0 _ 1 — e Ej=r la i 1 i U 0,l+Ei=i^l a i,i'"l,0 + 2 a i.i( M 0,l+ U l 1 o) ) i = 1,2, ...,N, (1.3) 



which evidently approximates (1.1) and (1.2). Here v? = u 3 (n,m), j = 1,2, ...N, is a set 
of the field variables. The subindex indicates shifts of the arguments n,m as follows u\ k := 
v?{jx + i, m + k) . 

The main result of the present article is in the formulating the conjecture below and proving 
it for numerous examples of Lie algebras. 



Conjecture, a) If A is the Cartan matrix of a semi-simple Lie algebra then (1.3) is Darboux 
integrable, in other words, it admits a complete set of integrals in both directions. Roughly 
speaking Darboux integrable difference-difference equations can be reduced to ordinary difference 
equations. 



b) If A is the Cartan matrix of an affine Lie algebra then (1.3) is S'-integrable (it admits exact 
multi-soliton solutions, generalized symmetries, Lax pairs, etc.). 

Recall that a function F (function G) is called an m-integral (respectively an n-integral) if the 
following identity holds D m F = F (or, D n G = G), where the shift operators D m , D n act due to 
the rules D m f(n,m) = f{n,m + 1), D n f{n,m) = f(n + l,m). Equations admitting a complete 
set of integrals in both directions are called Darboux integrable. 

Part a) of the conjecture is proved for the cases: An, Bn, Cn, G2, D 3 by finding complete sets 



of integrals. In the cases A2, B2, G2 the exponential system (1.3) reads as follows 



<1 -Li.l ^-,.1 ^ 0,1 ^,,1 _'>,2 



e - u i,i+ u i,o+ u o,i _M o,o _ x = e u o,i + "i,o _ ' u i,o ) (14) 

e -"i,i+«i,o+ u o,i _w o,o _ 1 = e - cw o,i +M o,i+ u i,o 

where c takes only three values c = 1,2,3. Here the values c = 1, c — 2, c = 3 correspond to 
the algebras A 2 , B 2 , G2 respectively. Let us introduce a notation for the nonlinear second order 
difference operator by setting 

A(n) := e -«M+ M i,o+«o,i-«o,o _ i j (1.5) 



then the last system takes a compact form 



A(u 2 ) 



g'^O.l+^l.O u l,0 



;i.e) 



Integrals in both directions for system (1.6) are given in Section 2 



Since the map defined in (1.3) converts any N x N matrix to a system of difference-difference 



equations one can easily specify the form of systems corresponding to any Cartan matrix canoni- 
cally related to a finite or affine Lie algebra both of classical series or exceptional. For instance, 
for the algebra -D3 with the Cartan matrix 



/ 



.4 



V 



we have an integrable system of the form 

AC* 1 ) = 

A(u 3 ) 



-1 



2 






2 



g — "o.i+^o.i+^i.o 



x o,i" 1 " u o,i-'" li i,o 



(1.7) 



Complete sets of integrals for this system is given in Section 2. 
For the Kac-Moody algebra A 2 with the Cartan matrix 



the corresponding system is 



A 



Aiu 1 ) 
A(u 2 ) 



2 

-4 



g u 0, !+■"!, u l,0 



e - 4 «o,i + "o,i+ M ?,o. 



;i.9) 



1.10) 



For the difference-difference systems corresponding the algebras A 2 , A[ \ A^ generalized 
symmetries are found, for A 2 the symmetries are local, for A± \ A^ they are nonlocal. In the 
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literature generalized symmetries are regarded as a criterion of integrability. Thus one concludes 



that systems (1.10), (3.1) provide new examples of 2 x 2 S-integrable quad graph models. 



Stress that the number of the field variables in all three systems (1.1), (1.2), (1.3) coincides 
with the rank of the corresponding Lie algebra. 

It is known that any system of difference-difference equations of hyperbolic type admits a pair 
of characteristic Lie algebras which are effectively evaluated (see [21] )• Characteristic Lie algebras 



of system (1.3) are closely connected with the Lie algebra canonically related to the given Cartan 
matrix A. In our opinion this relation could allow one to answer the second question in [TJ] (see 
the list of questions above). 

The problem of developing discrete versions of the Drinfel'd-Sokolov formalism is also chal- 
lenging since today systems of discrete equations are very popular. They have a large variety of 
application in theoretical physics, in discrete geometry, in the theory of tau functions in lattice 
Toda field equations. The present article would also provide an "experimental" background for 
creating a discrete theory parallel to the Drinfel'd-Sokolov formalism. 

The article is organized as follows. In Section 2 a method of discretization of Darboux in- 
tegrable models preserving integrability is discussed. The method is explained in details by the 



example of the Liouville equation. Then it is applied to the system (1.2) to get its difference- 
difference analog. The case B 2 is studied in more details. These simulations allowed us to guess 



the formula (1.3). Complete sets of integrals for the systems A 2 , B 2 , C 2 , G 2 , D 3 are presented. 

Generalized symmetries are evaluated for the systems corresponding A 2 , A±, A 2 in Section 
3. In Section 4 the concept of the characteristic Lie algebra (see [21]) of discrete models is briefly 
discussed. Here a complete description of such algebra is obtained for the system of difference- 
difference equations A 2 . 

In the fifth Section cutting off conditions for the Hirota equation preserving integrability are 
studied by using the method suggested in [2T] . It turned out that some of the systems in the class 



(1.3) can be obtained from the Hirota equation by imposing proper boundary conditions. Since 
the cutting off conditions are compatible with the Lax representation one can derive together 
with the reduced discrete system also its Lax pair. In §5 the Lax pairs for the difference- difference 
systems An, Bn, Cn, A^\ Dffl are found. 
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In §6 an algorithm is suggested to look for integrals via the Lax pair. In §7 widely known 
periodical reduction of the Hirota equation is discussed. The Lax pair for this system found in 
is rewritten in terms of the Cartan-Weyl basis for the algebra A$. Emphasize that the 



difference-difference system obtained as a periodical reduction and that given by the formula (1.3) 
corresponding to are not equivalent (see Remark 3 below) . 



2 Systematic approach to the problem of discretization of 
the Darboux integrable systems 

The problem of finding integrable discretizations of the integrable partial differential equation is 
very complicated and not enough studied. The same is true for evaluating the continuum limit 
for discrete models. In [22J an effective algorithm of discretization (as well as evaluation of the 
continuum limit) of Darboux integrable equations is suggested based on the integrals. In this 
section we discuss the essence of the algorithm and apply it to exponential type systems. 

2.1 Explanation of the method with example of Liouville equation 

Its well-known that the famous Liuoville equation u xy = e u admits integrals in both directions: 
= u xx — 0.5u 2 and W = u yy — 0.5w^. Indeed, it is very easy to check that D y W = and 
D X W = 0. 

Consider now the problem of finding all chains of the form 

t{n + 1, x) — f(x, t(n, x), t(n + 1, x), t x (n, x)) , (2.1) 
having the function I = t xx — \t x 2 as their n-integral. Equality D n I = I implies 

fx + ftt x + /ti/ + ft x t xx — -f = t xx — -t x . (2.2) 



By comparing the coefficients before t xx in (2.2) we have f tx = 1. Therefore 



f{x,t,t x ,t x ) = t x + d{x,t,t 1 ). (2.3) 



We substitute (2.3) into (2.2) and get d x + d t t x + d tl t x + d tl d — \t x 2 — dt x — \d 2 = —\t x 2 . 



or equivalently, d t + d tl — d = and d x + d tl d — ^d = 0. We solve the last two equations 
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simultaneously and find that d = e tl K(x,t\ — t), where K = Ce a^ 1 *) and C is an arbitrary 



constant. Therefore, chain (2.1) with n-integral I = t xx — \t x 2 becomes t\ x = t x + Ce^ tl+t ^ 2 . The 



last equation in its turn admits also an x-integral F = e^ 1 + e^ 1 * 2 )/ 2 . 
The next step consists in describing equations of the form 

v(n + 1, m + 1) = f{v(n, m),v(n + 1, m), v(n, m + 1)) 



(2.4) 



with m-integral I = e Vl v + e l 



Denote w 9 - 



j = 0, 1,2, and W\ := e U1 . In the new 



variables I = w ^ 2 is an m-integral of equation w^i = g(w,Wi,Wi). D m I = I implies 



iv 2 + w gi + wi 



Wi 
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(2.5) 



We differentiate both sides of (2.5) with respect to w 2 and apply the shift operator -D n 1 , we have 



1 



9i 



d: 1 ( -) = d? 

WiJ 



9i 



W-2 



9wi 



Wi 



w 



Therefore, 



WxW% 



w 



9 W 2 \ 

+ c{w,w 1 ), g x = \-c{wi,g). 

Wi 



We substitute (2.6) into (2.5) and get 



w 



c(wi,g) +w 1 . 



(2.6) 



(2.7) 



Substitution of (2.6) into (2.7) implies that c(w,Wi)w = c(wi,g)wx, or the same, c(w,Wi)w 



D n (c(w,wi)w). Suppose that equation w^i = g(w,wi,wi) does not admit an m-integral of the 
first order, then c(w,Wi)w = D n (c(w, wi)w) = C = const. Thus, c(w,Wi) = C/w. Finally, 

Cw^ 1 . Therefore, the equations (2.4) with m-integral J = e Vl ~ v + e Vl ~ V2 



g\w,w 1 ,w li 



becomes e~ Vl ' 1 ~ v = C + e~ Vl ~ Vl , where C is an arbitrary constant. Note that this equation is 
symmetric with respect to variables v± and V\. Therefore, n-integral for the equation can be 
obtained by simply changing in m-integral variables Vj into variables Vj, j — 1, 2. 
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2.2 Application of the algorithm of discretization to the systems of 
series B2 



Apply the reasonings above to the system (1.1) corresponding to the algebra B 2 . The first step 



has already been done in [27], where the differential-difference system was found: 



„.2 _ 2 _ -2« 1 +n 2 + U 2 
u l,x a x ~ c ■ 



(21 



admitting two independent n-integrals having the form: 



h = 2 «L + vL ~ 2 K) 2 + 1u\u\ - (u 2 xx ) 2 , 

+ulx(yL - uL) + u lxul(ul - 2u 2 x ) + «) 4 + {ulfiulf - 2(i4) 2 i£ 



(2.9) 



It was shown in [27J that system (2.8) admits also x-integrals: 

p 1 — e -«o+«i _|_ e - u \+ u h+ u 2- u l _|_ e u i-«2-«?+ u i _|_ e u \- u l^ 

The goal of this section is to construct a difference- difference system having the same n-integrals 
Fa and F 2 . 

To make the formulas shorter we change the variables: 

a = e~ u \ b = e~ u2 (2.10) 



where a, b are new unknowns. 

The given n-integrals in these variables read: 

a ai& 3 , a 2 6i a 3 

^1 = 1 7- H r- H , 

at a 2 b 2 a x b 2 a 2 

b a 2 b 2 a a 2 a 2 b x 63 ' 

t2 ~ T 2iT + 2 — 2" 2T r 7-. 

Oi ai 01 a\ z afb 2 b 2 

Here and below we omit the second index, because its value is zero for all considered variables. 

Now we will look for the equations desired in the variables a, b in the form: 

0-1,1 — f{ a o,o, o\ o, Oo,i, bo,o, b\ o, bo,i), 

(2.12) 

°i,i — 9( a o,o, a i,0) a o,i' ^o,0) ^i,0) &o,i)- 



Substitute the integrals (2.11) into the equations for n- integrals: D m F\ = F x , D m F 2 = F 2 and 



bring them to the following form: 

Qq,i fDjg gDnf Dlf_ = «o (hh «2&i <H 

f D n fD n g fD n g D n f a x a 2 b 2 a x b 2 a 2 

6q,i , ao,i 2 Dn9 9 ^iAJ ffAJ 2 _ 60 Qo 2 6 2 9 «o«2 «2 2 6i &3 
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f 2 9 



P 



f 2 D n g D n g b x a x 2 bi 



a\b 2 



By differentiating equation (2.13b) with respect to the variables 03,63 we get equations 

dg 



-Dl 99 



D n g oai aa\ 

— p> 2 — 1 fc» dg = 1)0,1 
D n g n dh b 2 dh b ' 



which imply immediately: 



g(ao,o, Oi,0) 00,1, 6o,q, 61,0, 60,1) 



^1,0^0,1 . <?l( a 0,0, a 0,l> ^0,0) ^0,l) 



Similarly, differentiation of (2.13a) with respect to 03,63 yields 

1 



:D 



2 df 



D n f n dai a 2 dai ao 



fD n g 



+ 



di df 



1 



b 2 D n fD n g D n f n dbi a 2 b 2 db x a b 
From these equations we get: 



D n\ a of ~ aiao,i) 



t( u u u \ — a i,o a o,i Cb\$ 

J\flO,Qi Ol,0) a o,i, °0,Q, 0l,0) °Q,lJ — 1 • 



a 0,0 



a 0,0 



By differentiating (2.13) with respect to a 2 we obtain 

ai,i 9 D n g x 



61, 



61 



2 Qq,i 2 02,161,1 1 



6262,1 <9a 2 02,1 0162,1 ai6 2 
2 



5 D 2 n _ a a 2 h 
01,101 0101,162,1 62,162,0 Ca 2 01 af0 2 



The last two equations can be rewritten as follows: 

d D n g 1 g x 



dct\ CL\ 1 O()0o 1 

— D n g x = 2 . 

aai Oo,iOq 



(2.13a) 
(2.13b) 
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By getting rid of the variable d ai D n gi in these equations we find 

D n g\ 9i 



°i,i a o,i 



^ 9i — Ci a o,i- 



The last implication is due to the assumption that F\ and F 2 are integrals of the lowest order. 
As a result we get equations: 

a i,i a o,o ~ a i,o a o,i = Cbifi, 

&1,1&0,0 — ^1,0^0,1 — Ci a o,i - 
After some rescaling we can put C = C\ = 1. Finally, we end up with the system of discrete 
equations 

a l, 1^0,0 — Ol,0 a 0,l = Ol,0j 
^1,1^0,0 — ^1,0^0,1 = i- 



corresponding to the algebra B 2 . Turn back to the original variables m 1 ,^ 2 (see (2.10)) 



Au l = e u o,i+ u i,o—"i,o 



Au 2 = e - 2u o,i+ u o,i+«?,o. (2.14) 



By construction the last system admits pair of integrals (2.9) in which the second index for the 
variables w 1 ,?/ 2 is omited, since its value is the same for variables. To find the integrals in the 
other direction one uses the discrete symmetry n — > — m, m — > —n of the system. 

2.3 Integrals of the systems corresponding to the series A2, G2, 
2.3.1 System corresponding to the algebra A 2 

The Cartan matrix of A 2 is: 

A = 



2 -1 
-1 2 



The discrete system for A 2 looks as follows: 



Aiu 1 ) = e "o,i+"i,o-«f >0) (2.15) 
A(u 2 ) = e~ u o,i +u o,i+ u ?,o. 



11 



For system (2.15) m-integrals are 



P-y — g M 0,0+ n l,0 _|_ g U 0,0+ U 1,0+ U 1,0 M 2,0 _|_ g U l,0 n 2,0 
— g~ M 0,0+ M l,0 -)- g n 1.0 — n 2,0~ U l,0+ U 2,0 _|_ g u 2,0 — "1,0 _ 



2.3.2 System corresponding to the algebra G 2 

The Cartan matrix of G 2 is: 

A = 




The discrete system for G 2 is of the form: 

A(m 1 ) = e u °> 1+ "i.°~<° 
A(w 2 ) = e _3M o. 1+?t o. 1+ "i' . 

Its m-integrals: 

P^ = g u l,0 — u 2,0 _)_ g u -l,0~ U -2,0 _)_ g u l,0+ u l,0 — M l,0 — u 2,0 _)_ g n -l,0+ n 0,0~ u 0,0~ "-1,0 

_|_g u i,o+ 2 ' !1 o,o~ u o,o — 2u i,o _j_ g u o,o+ 2tl o,o~ u i,o — 2u -i,o _|_ 2e 2u o,o — U -l,0~ U l,0 ^ 

^ = g u 2,0 _u 3,0 _|_ g 3n 0,0+ n 0,0 _u l,0 _3u -l,0 _|_ 3g u 0,0+ tt l,0 _u 2,0 _u 0,0 _)_ 3g u 0,0+ 11 l,O~ u -l,O _,1 2,O 
_|_3g u ?,0+ 3?1 M~ 2u 0,0- u 2,0-"2,0 -|- 3g 3u l,O+ U 0,0- 2u 2,0 _|_ 3g 3u 0,0+ n J,0-' u -l,0 

-j-3g 3u o,o+ u i,o _u o,o~ 2u i,o _ "-i,o _|_ 3g u i,o+ u i,o — "a.o - u -i,o _|_ 2e 2u i' 0+ "o,o~ u 2,o 

_)_g 3u i,0+ n 2,0 _tt l,0 _3u 2,0 _)_ g 3M l,0+ 2 "l,0 _3n 0,0 _2u 2,0 _)_ g 3u 0,0+ 2u l,0 _2u 0,0 _3?t l,0 _)_ g U 2,0+ U -l,0 

2.3.3 System corresponding to the algebra .D 3 

The Cartan matrix of D 3 is: 



(2.16) 



(2.17) 



A 



( 



\ 



2 -1 


-1 1 


-1 2 





-1 


2 > 



The discrete system for D3 looks as: 



A(u 2 ) = f 
A(u 3 ) 



g u O,l+ M l,0 u 0,l n o,i ; 



(2.18) 



g«0,l+ U l,0 — U l,0_ 
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Its m-integrals: 

= g U 0,0 — u -l,0 _|_ g u l,0 — u 2,0 _|_ g u l,0+ u l,0 — u 2,0~ n 0,0 _|_ g u l,0+ u O,0 — u 0,0 — u l,0 
_P 2 = g u O,0~"-l,0 _|_ g u l,0 — u 2,0 _)_ g"l,0+"l,0~ u 2,0~"' u i _|_ g U l,0+ M 0,0~ ^.O - u l,0 

= g"l,0 — u 2,0 -j- g u O,0 — u -l,0 _)_ g u O,0 — U l,0+ U 0,0 n -l,0 -j- g u O,0 — M -l,0+ M 0,0 — u l,0 
_)_g M i,o — "o.o+^o.o - M l,0+ n O,0 — U \,Q _|_ g U 0,0 — U l, O+^O, — U -l,0+ U 0,0 — u — 1,0 _ 

To find the integrals in the other direction one uses the discrete symmetry n — > —m, m — > —n 
of the system. 



3 Symmetries of discrete systems 

In this section we demonstrate that the discrete systems admit generalized symmetries. For the 
case of the simple Lie algebras the systems admit local symmetries while for the case of affine 
algebras the symmetries are nonlocal. 

3.1 Generalized symmetries for the system 

Higher symmetries for system corresponding to A 2 are found by using the method suggested in 
[29] . (see also (SO]). The lowest order symmetry is of the form: 

du o,o _ 3e"o,°""i.° _ y 



dt " D^Fi 



du?. 



0,0 - I _ 3e"0," 



2 2 

Wn n — u_- 



dt - D~ 1 F 1 ■ 

The next order symmetry depends on arbitrary functions f(x,y), g(x) of two and respectively one 
variable and has the form: 



/(Fl,^)+ g (^F0 e < -uL 1|0 + D -lp 2) _ g(Fi) + g( D -l Fl ) 



du 0,0 

dt D^F 2 

f(F 1 ,D^ 1 F 2 )+f(D- 2 F 1 ,D- 2 F 2 )-g(F 1 ) u \ -u\ „ 

du lo _ f(FuD^ 1 F 2 )+f(D^ 2 F u D- 2 F 2 )-g(F 1 ) u \ Q -v?_ x „ _ f( r>-2 p r>-2 p \ , 
~~" o; 1 ?! ^ \ m 1' m 2Jt 

f{D m 1 F 1 ,D^ l 2 F 2 )+g(D- 2 F 1 ) u 2 -u\ „ 
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3.2 Evaluation of a generalized symmetry for the system 

The method for searching nonlocal symmetries is illustrated with the following example: 

g— «1,X— "0,0 g— «1,0— "0,1 _ g— 21)1,0 



g— Dl,l— D0,0 g— Dl,0 — D0,1 _ g— 2«0,1 



corresponding to the algebra with the Cartan matrix: 

2 -2 
-2 2 

Recall that the continuous version of the system is: 



(3.1) 



A 



u 



exp(2w — 2v), v xy = exp(2t> — 2u). 



(3.2) 



It can be proved that system (3.2) does not have any local generalized symmetry, however it has 



nonlocal generalized symmetries, the simplest one can be represented in terms of s — u + v, r 
u — v as follows: 



- 2r 3 x , s xt = 2r x r xxx - 3r x - r 2 xx + F(W, W x ) 



(3.3) 



Here F is an arbitrary function and W = s xx — r 2 is the ^/-integral of system (3.2). Note that 
years ago in |9] it was observed that equation r f = r xxx — 2r x is consistent with the system (3.2). 



For the fully discrete analogue (3.1) of the system (3.2) we have a very similar situation 



The system (3.1) does not have any local generalized symmetry. One has to look for a nonlocal 



symmetry. The main trouble arising in this case is connected with the guessing of the structure 
of non-locality and the form of the symmetry. 



Rewrite system (3.1) in the form: 



a i,i a o,o ~ a i,o a o,i — bi , 

&1,1&Q,0 — ^1,0^0,1 — ^0,1' 

where a id = exp(-u iii ), b id = exp(-u i)i ). 



(3.4) 



We will look for the generalized symmetry of (3.4) in such a nonlocal form: 



dtCllfl — /l(O0,0) fll,o)^i a O,0 + /(O-l,0) &0,Cb a l,Cb ^2,0; ^-1,0) &0,0> &1,0> &2,o)j 
dt^lfi = PlXVo) ^l,o)<9^0,0 + <?(°-l,0> a 0,Cb a l,0) a 2,0, &-l,0) ^0,0; ^1,0; &2,o)- 



(3.5) 
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Here fi,gi,f,g are unknown functions. 
From (3.5) we can find: 



d t a 1A = f m (/i)9 t a ,i + D m f, d t b 1A = D m (gi)d t b 0tl + D m g. 



(3.6) 



After differentiation (3.1) with respect to t by means of the systems (3.5), (3.6) we get: 



(3.7) 



•9^0,0(01,1 - a ,i/i) + 9 t a ,i(ao,o-D m /i - «i,o) - 2<9i& , fri,ofi , i = • • • , 
-2a ,i<9 t a ,i + ^60,0(61,1 - b 0>1 gi) + d t b 0>1 (6 o ,oAn0i - &i,o) = • • • 
Here the right hand sides do not depend on the derivatives of dynamical variables with respect to 



t. By applying the operator D n to both sides of (3.7) we obtain: 



d t a 0fl fiD n (a lt i - a ,i/i) + d t ao tl D m fiD n (a 0)0 D m fi - ai,o) - 1dtb\$b 2 $g\D n g x 
-2ai + d t b^g x D n {b x;v - 6 o ,i0i) + <9^o,iAnfi , i£ , n(&o,oAr t 5'i - 61,0) 



(3i 



System of equations (3.7), (3.8) is a system of linear algebraic equations with unknowns d t aot 



dtbo,o, d t a 0: i, d t b 0) i. If the determinant of this system is different from zero, then due to the 
Cramer's rule the system has unique solution, and therefore the searched generalized symmetry 
is local. This output is in contradiction with our previous study proving the absence of local 
symmetries. Thus the determinant should be zero: 

— fl a 0,l O0,0-Dm/l — a l,0 — 26i,o(7l 

/i(a2,i - a>i,\D n fx) (ai >0 D n D m fx- a 2 ,o)D m f 1 -2b 2>0 giD n g 1 

-1a>\,\V m h giih,! - &i,i-Dn0i) D m g 1 {b Xfi D n D m g 1 -b 2 fl) 

-2a ,i 61,1 - 5-160,1 b ,oD m gi - 6 1)0 

(3.9) 

Vanishing of the determinant implies 



A( a o,0) a i,o) — a i,o/ a o,0) gi(bo,o, bifi) — b 1)0 /bofi, 



(3.10) 
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In virtue of (3.10) equations (3.7), (3.8) take the form: 



dtdo,o . d t a i <9t&o,o 2 ao i , a^o „ „ 
+ 2 ^ — = r~ 9 + 7^ ~ ut D ™L 

00,0 °0,1 °0,0 01,0 °l,0 °1,0 

d t a 0fl d t a 01 O 5 t 6 , 2 ai , a 1Q / D m f 2 

H 2- = - — D n g + —D n f - ^DnDmf h j—g, 

0i,i 0i,o 



00,0 



ao,i 



0o,o 



02.0 02 



2.0 



0-1,0 



5 t 6 ,o , <9<0o,i o «9i0O,l 0o,i , 0o,o n 

+ 2 = -=-p + -^D m g, 



(3.11) 



00,0 00,1 00,1 

d t b ,o . 5*00,1 5t0o,i 0i,i n . 0i,o n n 2 
+ 2 = —D n g + —D n D m g H 

T i,l 



A o,i 



l o,i 



0o,o 0o,i °o,i a i,i a i,i 0i,i 0i,o 0i,i 

We can see that the left hand sides of the first and third equations as well as the left hand sides of 

the second and fourth equations of the last system are identical. Comparison of their right hand 

sides gives a system of two equation for /, g: 



0i,o n n , 0i,i n p 0i,o0o,i n - 0i,i0o,o n 

— g L) n L) m J --j J-JnJ 7~2 ^rnj + — ~pi J 

°2,0 2,0 01,101,0 



01,001,0 



^2,0 



-D n g = 0, 



0i,o n n . 0i,i n 0i,o0o,i n 0i,i0o,o ,2 
-^-D n D m g + -^-D n g - ^£> m p - ^g H D m / = 



n,i 



"i,i 



0i,i0o,i 



0i,o0o,i ' 



0i,i 



Differentiate (]3.12|) and (]3.13|) with respect to 03,0, 63,0, 6-1,0 or a_i 5 o and get: 
0o,o0o,o df 



0?o 56_ lj0 



OO,O0Q,O 5/ 

96-i,c 



A o,o 



5/ 



n.o 5/ 2a_i,ia 01 



a i,o 5a~i,o / a l,o 5a-i,o a l,l°0, 



D r 



(3.12) 
(3.13) 

5/ 



01,0 dg _ ai,i 26i,o6 2 ,o n 5# &i, &o,i 5# 

0o,o 5a 2 ,o 0o,i 5a 2 , 



(3.14) 
(3.15) 



df 26i )0 5# _ 01,100,0 5/ 



5a2,o 00,1 <9a 2 ,o 0,1,000,1 m 5a 2 ,o' 56 2 ,o a ,i 56 2 



56 2 , ' 01,001,1 m da 2 ,o 0o,o0i,i 56 2 ,o : 
df 2b lfi dg _ 2a 0fi b 2 fi df Qq,qOi,i 5/ 



01,000,1 



9 5 | 26 1|0a § [0 g/ = g ! og ! -i £) - 1 ^ 



56-i,o 0i,o0i,-i0o,o 56_ 



1,0 



2 w^r m 50-1,0 



5a 2 , ' 00,101,0 m 562,o' 
(3.16) 

(3.17) 



dg 2o§ [0 6i [0 df _ 6j i0 o _i& _i x 4a Q i i0 a_ li o6 li o fl/ 2o_i t o 5# 



5a-i,o 01,061,-160,0 5a_i i0 00,060,061,- 



5a_i )Q ai,o6 -161,-160,0 56-i,o &o,-i 56_ lj0 ' 

(3.18) 

We use the last system of eight equations in order to specify dependence of the functions /, g on 
the variables 0-1,0, 02,0, 6-1,0, 62,0 corresponding to the lowest and highest values of the first index. 



From the first equation in (3.14) we get: 



df C 7 a? 



56_i ; o 6o 5 o0o,o 
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(3.19) 



Then the second equation takes the form: 



D r 



A 0,0 



Qq,o df 



df 2C 7 a_i 5 ia 0i 



^0,1^0,0 



(3.20) 



It can easily be proved that the expression 2a bo \ll° Q 1 does not belong to the image of the operator 
D m — 1, so we have to put CV = and then 

df c 4 4o 



From the pair of equations (3.15) we have: 

dg 



(3.21) 



da 2 n 



0. 



dg _ C 2 b 0fi 

db 2 n 6i 



(3.22) 



Continuing this way one can find dependence of the functions /, g on the variables a_i,o, 02,0, 0-i,o> ^2,o : 
1,0^2,0^ . n ( a ,o«2,o . a o,o&|,o 26 0j0 ai i0 62.o\ . a i,o a -i,o „ , , u u \ 

bi + Lj 2 I 1 "2 1 I H 2 ^4 + JnOo,0, Ol,05 °0,0, °l,0j, 

°1,0 \ a l,0 a l,0°l,0 °1,0 / a o,o 

^1,0^-1,0 , &i,o a -i,o . 26 lj0 a li0 a_i i0 \ fef j0 a_i )0 



a 0,l 

= ~~6 
bo,ob 2) o 
bin 



C 2 + Ca 



ft 2 

u 0,0 



+ 



fl 2 ft 2 

U 0,0 U 0,0 



+ 



l 0,0 



+ 



a 0,0^0,0 



-C5 + gi{ao,o, aifl, &o,0) °i,o) 



(3.23) 



For the functions f\ , g\ from the same system ( 3. 12 ) , ( 3. 13 ) we obtain the following non-homogeneous 
equations: 



Q i,o n n , n f 0-1,000,1 n , 01,100,0 , 2 

JJ n JJ m Jl — -g— iA/l To" i/m/1 H To - Jl — 7 -^n^l 

°2 01,1020 Oi,q0 1 02,0 



6 2 
"2,0 



D n D m gi + 1 -^D n g 1 



bi,obo, 



h,ib 0t o 



"1,1 



"1,1 



1 2 J -D m g 1 — — J — 2^gi H -D m /i 

°l,l a 0,l °l,0 a 0,l 



2 

01,1 



(3.24) 
(3.25) 



Differentiating (3.25) with respect to the variable a 2 , we get: 

2C 2 ( a o,l a l,0 — °1,0^0,1 



00,0 dgi ( a , dgi 

u 7 



60,0 da\fi \°o,o daifi 

From this equation we can find: 



Ca — Co 



Oo,oOo,lOl,lOl,o 

Qq,o dgi _ 26*201,060,0 
00,0 <9oi,o o ,o6i,o 



+ 2(C 4 - C 2 ) 



l 0,lO-l,0 — °1,1°0,0 

oo,oOo,oOo, 160,1 



+ c 6 , 



therefore 



91 = a 



a l,oK,0 , ^ Oi n6o,0 



a o,o & i, 



+ C 6 



a 0,0 
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+ £2(00,0, &o,oA,o)- 



Differentiating (3.25) with respect to 62,0 we get: 
dg 2 



therefore 



_ dg 
dbioj <%i,o 



2 2(C 1 -C 6 )^5. 

flo,o a o,i 



2(d - C 5 



a 0,l a 0,0 a -l,l — ^0,0^1,0^0,1 



a O,oVo a O,l&0, 



(3.26) 



C*6 — Cl> C5 — C\ 92 — £-3^1,0 + ^(^O) ^0,1 



Now apply the operator D n x to the equation (3.24) and then differentiate the obtained result with 
respect to the variable 6_i : 



D 1 ^0,061,0 dfi ^ _ 00,061,0 d/i _ 2(7 2 & i< 0&0 ' 1 ~ ai >° a o,i 



a i ,0/ Q i 



a o,oOo,i6i,o6 



1,1 



(3.27) 



therefore 



„3 Z,2 „2 z. 

/ _ n a i 5 o°o,o . ^ a i,o°o,o , , , s 

Jl — ^2 2 l2 1" ° 7 l r 721^0,0, ^1,0) °l,0j- 

a 0°1 a 0,0Ol,0 



(3.28) 



Substitute in (3.24) the expression found instead of fi and differentiate with respect to 02,0 to get: 

a i,o a o,i ~ ^1,0^0,1 



D r 



dai 1 



2(C 1 -C 7 )- 



Oo,o a o,i6i,o6 



1,1 



therefore: 



(3.29) 



CV — Ci, / 2 — C 8 ai,o + /3(ao,0) &i,o)- 



For the functions /3, #3 we have a system: 

ai >° n n / Q 1>1 n f 01,0^0,1 



ai,iao,o 



, » D n D m f 3 — T^D n fz ' D m f 3 H ' / 3 — - — D n g 3 — 2(C 8 — C 3 ) 

bin bi n ai.ifes " '- 



y 2,0 



2,0 

&1,0 n n &1,1 n &1,0&0,1 n 

-^—jy n jy m g 3 -\ — ^-JJ n g 3 — ^—u m g 3 



a l,o^i,o 
61,160,0 



2 

62,0 



a 



1.1 



a 



1.1 



6i,i a o 1 



g 3 + D m f 3 = 2(C 8 -C 3 ). 



6i,o a o 1 °i,i 



We see that (3.30) is satisfied identically if one chooses 



(3.30) 



C 8 = C 3 , f 3 = 0, g 3 = 0. 
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Thus we find the final form of the symmetry searched: 
dtCLi,o _ d t a Qt0 / ao,o^2,o a i,o^o,o 



a 1.0 



a o,o 



°1,0&1,0 a 0,0^1,0 

+ 



/ 2 t.2 2 t.2 

/ 0x^0-1,0 26o,0^2,0 a 1.0"0,0 ^0,0^2,0 a 0,0"2.0 

+ ^2 | 2 ' ^ ' " — ^ ' R ' s — 77 "~~ 



l 0,0 



/r 

°1 ' 



n 2 h 2 
a 0,0°l, 



'1,0 



« 2 /) 2 



3- 



dtbtfl dfbop / 0-1,0^1,0 _|_ ai,o^o,o 



&i,o 



&0,0 



(3.31) 



+C 2 



a 0,0^0,0 a 0,0^1,0 

2ai )0 a-i,o &o,o^2,o . a i,o^o,o . &-i,a&i,o . a -i,o&i,o 



+ 



+ 



a 6,o ^i,o a o,o^i,o 



+ 



b 2 
°o,o 



+ 



n 2 h 2 



C 3 



Here Ci, C2, C3 are arbitrary constants. 

By analogy with the continuous case we have a differential constraint completely consistent 



with system (3.4), which is obtained from the last system by applying the operator (D n — 1) 1 to 



the difference of the two equations: 



bo,o 0-0,0 a o,obo,o 
Under the Cole-Hopf type transformation 

«o,o = 



a -l,0 a l,0 , &-l,0&l,0 . a -i,o^i,o 



l o,o 



h 2 
°o,o 



a 2 h 2 



«0,0 



bos 



&0,0 



(3.32) 



(3.33) 



system (3.4) converts to the following one: 

^1,100,000,1 — Oi 5 oao,i^o,o = &i,o(^o,o — Oo,i)j 

^1,1^0,0^0,1 — &i,o^o,i^o,o = ^0,1(^0,1 — ^0,0) • 
This transformation brings our symmetry to a local symmetry of the new system (3.34): 



(3.34) 



<9 t a ,o = Ci &i,o + 



0,0 



4),o 



C 2 



0, 



0,0 



+ 



2&i,oOo. 



b 2 
°o,o 

a-1,0 



+ a ,o + C 2 



a -l,0 

2o-o,o^o,o 
a-1,0 



+ 



a: 



0.0 



J o,o 



h 2 
°o,o 



+ 6i,o + f h 



ai,o 



6? 







ft 2 

°0,0 



^0,0 &-l,0 



+ 



°0,0 

bio 



+ C 3 a . 



0) 



(3.35) 



-1,0 



(2) 

3.3 Generalized symmetries for the system A 2 



The Cartan matrix of is: 



.4 
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The discrete system for looks as: 

A^ 1 ) = e <i +u l°- u l°, (3.36) 
A(u 2 ) = e _4u o> 1+u o. 1+u i.°. 

Generalized symmetry 

(yi Q — ul ) t = Ci ^e 2u ifi- u lfi- u 2fi _)_ e 2u i,o- 2u o,o+ u i,o- u l,o _)_ e 2u o,o- 2u i,o+ u i,o- u o,o _)_ e 2u o,o-"i,o-"-i,o j 

+ C 2 ^3e 4U l.° -2U 0'° -2U 2,0 + e U l,0- u 0,0+ u 2,0- u 3,0 -|- 3 e - u 2,0+ u l,0+ u 0,0- u 0,0 

_|_3 e 4 "i,0- u 2,0- 3u 0,0+ u l,0-"2,0 _|_ 2e M 1.0- U 2,0+ u 0,O-"-l,O _|_ e 4u l,0- 4u 0,0+ 2u l,0- 2u 2,0 
_|_2e U 2,0 _2u l,0- u 0,0 _|_ 3gW},0- u 2,0+ u l,0- u -l,0 -|- e - u i,0+ u -l,0+ u 0,0- u -2,0 
+ e - 4 "i,0+ 2M io- 2u O,()+ 4n O,0 _)_ e 4u l,0- u 0,0+ u 2,0- 3,t 2,0-' lt l,0 -f e -"l,0+ U -l,0+ U 0,0- U -l,0 
_|_3g- 2 «l,0- 2 "-l,0+ 4u 0,0 _|_ e -«0,0+ u L)+ u 2,0- M 3,0 _|_ 3 e 3u l,0+ M io- u -l,0- u O,0+ 4?t O,0 

_|_ e -"i,o- u i,o- 3u -i,o+ u o,o+ 4u o,o^ _|_ (7 3) 

( M 2 o _ MQ ) t = 2Ci ^e _2 "i- o+ "i> ()_?t "' 0+2u 0' + 2e -u i>° - '"- 1 ' 0+2u 0' + e _u i. 0+2u O' 0+u o- 0_2u - 1 ' ^ 

+2C 2 ^2e u i.° _u 2,o+ u o,o- u -i,o + 2e _u 2,o +u i.o +M o." _M o.o + 4 e ~ 3u M+"L)- u -i,o- n o,o+ 4 "o,o 
4- e - 4u l,o+ 2u ?,o- 2u o,o+ 4u o,o + 2e u i.o _u 2.o +u i.o~'"- 1 .o + e ~ M 2.o +2 "i>«~ M »>" 
_l_ 2e _u i' 0+ulL:L > 0+?t(1, ' ~ M - 2 ' + 2e — "i.o+^-i.o+^o.o - u -i,o _|_ 4e — u i>o — u i,o — 3u -i,o+ u o,o+ 4u o,o 

+ 2e _ "i- 0+ " ll) > ()+?t "' 0_M - 2 ' + e _u i' 0+2u °' 0_u - 1 - + Ge -2 "!. 0-2 "- 1 - " 1 " 4 "^ 
+ e - 2 «?,o+ 4u i,o+ 2u o,o- 4u -i,o j + 2C 3 . 

(3.37) 

Here Ci, C 2 , C3 are arbitrary constants. And for combination we have a local constraint: 
(2uJ >0 - ul ) t = 2d (e 2 <°- u l°- u -^ + e H°- 2u -i.° + <°- u i.°) + 

2C 2 "-i.O - U 2,0+ U 0,0+ U l,0 _|_ g u l,0+ U l,0~ M -l,0~ n 2,0 _)_ g _3M l,0+ n l,0+ 4u 0,0~ U -l,0~ M 0,0 

+ 3e~ 2M i.o +4 "< L »>o~ 2M -i,o + 3e _u lo- u i,o+ 4u o,o- 3u -i,o+ u o,o _|_ e - u i,o+ u -i,o+ u o,o- u -i,o 

4- e - u i,o+ u o,o+ u -i,o- u -2,o 4- e - 2u i,o+ 4u ( 1 ),o- 4u -i,o+ 2M S,o _|_ e - u ?,o+ 2u o,o- M -i,o 

_|_ 2g~ U l> ^ U(1 )."^ M( W' — U -2,0^ 

(3.38) 
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4 Characteristic m— algebra for the case 

Let us describe briefly the properties of the characteristic Lie algebras of the system 



ao,o«i,i = ai,o«o,i + h,o, (4-1) 

0O,O°1,1 = °1,0°0,1 + a o,i- 



System (4.1) corresponds to the simple Lie algebra A 2 . First we introduce the notion of the 



characteristic m— algebra for the system (4.1 ). Recall that according the definition an m— integral 
F(ao t o, bop, a\ t o, b^o, a_i j0 , &-i,o> • • • ) should satisfy the equation D m F = F. In the coordinate 
representation this condition reads 

) = F{a ). (4.2) 



Evidently the right hand side in (4.2) does not depend on the variables a 0j i and &o,i hence the 



conditions hold ^D^F = 0, g^D^F = as well as Y X F = 0, Y 2 F = where Y 1 := 
D^Q^—D m , Y 2 := D^g^D m . Denote through L n the Lie algebra generated by the operators 
X\ = da ® , X 2 = db ® - , Ki, Y2- Here the first two operators are the first order linear 
differential operators, or vector fields. The last two operators can also be rewritten as vector fields 
of the form 

Y = - + ( SM. - ^A- 1 ^ 9 + ( a -i.o + fe o,-i ^ 19 
(9a 0) o \Qo,o a o,o^o,o a o,-i / <9&i,o \ <2o,o Oo,o a o,-i / «9 a -i,o 

19 / a _ 10 1 \ d 

+7 ~— - ^ + - + • • • , (4.3) 

00-1^01,0 \ a o,oV-i a o,oa ,-i/ o6_i j0 

cwo,o \ & o,o »o,o»o-i/ o6i, V & o,o 00,000-1/ o6_i j0 
Taking commutators of the vector fields X 1 ,X 2 ,Y"i, Y2 we get vector fields 

Pi = [X 1 ,Y 1 ],P 2 = [X 2 ,Y 1 },P 3 = [X 2 ,Y 2 ],P 4 = [X 1} Y 2 ]. 

Obviously the vector fields X%, X 2 , Y%, Y 2 , Pi, P 2 , P 3 , P 4 belong to the m— algebra. We use the 
following lemma to show that m— algebra is of finite dimension. 
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Lemma 4.1 Suppose that the vector field 



\ - I d d 

K = 2_^[ a k - h a- k 

k=i 



d 



(4.5) 



satisfies the equality D n KD n 1 - 
variables a and b, then K = 0. 



da fc oa_ k J f-f V <™fc c/o_ fe/ 

7 fc=i x ' 

ZiX, where h is a function depending on shifts and derivatives of 



One can easily check that 



D n Y x D~ 



1 _ a o,o 

11,0 ' 

1 oq,-i v _ ap,-i 



D n X x D~ v = ^Xi, D n X 2 D~ 1 = -^X 1 + ^X 2 , 



-y 2 , D n Y 2 D~ l = b ^Y 2 . 



Put Xi = a , Xi, X 2 = 60,0X2, Yi = ao-iYx, Y 2 = b -iY 2 , then 

D^D- 1 = Xi, D n X 2 D~ l = -^Xi + X 2 , 

n n ao,OOi,o ' 



D nYiD n ' V : fto.-V-i^ 2 ' 



DnYzD' 1 = Y 2 . 



Lemma 4.2 



D^D- 1 =P X 



"0,0 



bo, -161,-1 
61,0 



Y 2 , 



DnP-iD- 1 =P 2 + ^^P 1 - 

n z n z 00,001,0 



Q(),-l 

60,-161,- 



-P 3 + 



+ 



( «u,-i 61,0 _|_ 01,-161,0 \ x + ( a o,-i6o,o _|_ 01,-100,061,0 \ y 
\ Oo,o ai .o oo,oo? j0 / 1 V^o.-i 6 !.- 1 Ol > o6o,-l6f _ 1 J 2 ' 



DnPsD' 1 =P 3 



61,-1 
oo,oai,o 



Xl 



Lemma 4.3 For c = 1, 2, 3 we /lave P 4 = 0. 





x, 


x 2 




Y 2 


Pi 


P2 


P 3 










Pi 








Ri 





x 2 








P2 


P 3 


i*l 


R2 


-2X 2 


*i 


-Pi 












R 3 


R4 


F 2 





-P3 











R4 





Pi 





-Ri 


-2Yi 








Qi 


Q2 


P2 




-R 2 


-R3 




-Qi 





Q 3 


P3 





2X 2 


-R 4 





-Q2 


-Q 3 






(4.6) 
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_ 2oo,-i6g „ 2o 0l -i6g ^. 2o ,-i&o,o v 

Ji/; — — — j q — tt Jo T r v\ o i 

ao,0&o,-l ao,0&o,-i oo,o6o,-l ' 

_ 2a§ ,60,0 D 2a , -160,0 \> , 2a o,-i b O,o ^ 

ag b ,-l ao, 0»0,-i 1 ag &o,-i ' 

P 4 = -Yi ~ - ^P x , 

a o,o a o,o x ' 

Qi = - 3P 2 - P 3 + 4^Xx + 

oo,q6o,-i d a o,o°o,-i 60,-1 A ' 



Q 2 = g*Y 2 - 2P 



3- 



Q op 1 QQ,-1&0,0 "0,-1 Y a O,~l°0,0-y> 

Any m-integral P(. . . a 0j o, 60,0, a i,o, &i,o, ^2,0, &2,o • • • ) must satisfy equations 

X 1 (F) = 0, X 2 (P) = 0, Yi(F) = 0, F 2 (P) = 0, 

Pi(F) = 0, P 2 (P) = 0, P 3 (P) = 0. (4.7) 



Solving system (4.7), it is enough to assume that Pi depends on 6 ,o, 00,05 °i,05 fl i,05 ^2,05 02,0 



and P 2 depends on 0-1,0, ^0,05 fl o,05 01,05 a i,o ; °2,05 we find two independent m-integrals 

t-, ^0,0 . a 0,0^2,0 . a 2,0 T-, a -l,0 . a l,0^0,0 . ^2,0 , A Q x 

^1 = 7 1 7 1 5 -^2 = 1 7 h 7—- (4.SJ 

°1,0 «1,0°1,0 a l,0 a 0,0 a 0,0°l,0 °1,0 

5 Cutting off constraints for the Hirota equation and dis- 
crete Zakharov-Shabat systems 

It is well known that Hirota chain 

t 0,0 t 'l,l ( '1,0''0,1 ~~ 4,0 ''O.l 

admits the Lax pair consisting of two linear discrete equations 

,j _ r l,0 r 0,0 ,j jj _,j *0,1 %0 (52) 

ri,o — ^0,0 ^0,0 5 Wo,i — ^0,0 + j j ^0,0 • \°- z ) 

Vo n,o ^0,0*0,1 

Here the lower indices mean as previously shifts of the arguments, and the upper index enumer- 



ates the field variables P and eigenfunctions ip J . Exclude from the system of equations (5.2) all 



the eigenfunctions with the upper index different from j. As a result one gets a linear discrete 
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hyperbolic equation for ipi 



L l,l L 0,l , j , t 0,0 £ 'l,l 



yj+lyj 



^0,1 + ^TTT^o.o 

t l,O t '0,l 



0. 



Equation (5.1) is invariant under the transformation defined as n 



m, m 



this transformation the Lax pair (5.2) transforms to a Lax pair 



V-1,0 



2/o,o + 



-i,cro,o „j — 1 



-2/6,o 



% ,-i r o,o J 

*J+M 2/0,0 i/0,0 > 
'0,0 l o,-i 



J+l 



and equation (5.3) transforms to an equation 



t j t j+1 

f 'l,0 t '0,l j 

''Co fc i,i 



yj 

t 'l,0''0,l j 
t 0,0 , 'l,l 



yj yj+l 
t 'l,0 , '0,l j 

^0,0 l l,l 



0. 



(5.3) 

1 — n. Under 
(5.4) 

(5.5) 



Put yU 



-1,0 .J 



g 00 , then Lax pair (5.4) transforms to a Lax pair 



£-1,0 



j+i ^2/o,o + 2/o,o J > 



r 0, 0-2,0 



yo,-i 



C 0,-1 C -1,0 _7+l 



'O-l'-l.O J - - rr- ■ 

p+i^+i yo.o j+ifj+i 2/0,0 

r 0,0 1,— i 'o,o r -l,-l 



and equation (5.5) transforms to an equation 



t j t J+1 , , 

j J t, l,0 t -l,l j '1,0'-1,0 j _ n 

yi,i 2/i,o ,7+1,1 £0,1+ ,7+i,7 2/o,o — u - 



t j t j+1 



yi + lyi 

''0,1 ^o 



yj + lyi 
'UO t 0,0 



(5.6) 



(5.7) 



Thus we have two different Lax pairs for the Hirota chain and consequently two families (5.3), 



(5.7) of linear discrete equations enumerated by j. Study the question when an equation from one 



family can be related, by a linear transformation, to an equation from the other family. To this 
end we evaluate the Laplace invariants of these equations (see [20]). 

Recall that the Laplace invariants of a discrete hyperbolic type equation of the form: 



/l.i + Oo,o/i,o + &o,oio,i + Cq,o/( 



0.0 







(5. 



are given by (see 



fco,oCi,o 

^0,0^1,0 



fro,iCo,o 

^0,0 a l,0 



(5.9) 



By virtue of these formulas the invariants K\$,Kty and Ki g ,K 2g of equations (5.3), (5.7) are, 
respectively, 



K 



lib 



f J f yJ+ly. 
t 'l,0 t 'l,l £ '0,l L 



2,1 



K 



lg 



+3 +3 jj'+ljj+l ' 

t, o,o fc 2 ) i fc i,i n,i 

t 0,0''2,l , '0,0 fc 0,l 

jj jj jj'+l jj+1 ' 
t l,0''l,0''-l,0 t 'l,l 



K- 



2t 



K- 



2S 



yj yJ 
^1,0 6 0,1 

yJ fJ ' 

( '1,1''0,0 
fc 0,0 ''-l,! 
t -l,0 t 0,l 



(5.10) 
(5.11) 
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It is known that two linear hyperbolic type equations are related to one another by a linear 
transformation only if their corresponding Laplace invariants are equal. Evidently in generic case 
coincidence of the Laplace invariants generates two constraints on the field variables P = P(n, m). 
Only for some special cases it gives only one constraint. We are interested in such special cases. For 
instance pair of equations Ki^(n, m,j) = Ki g (n + 1, m, j — 1), K 2 ^(n, m,j) = K 2g (n + 1, m, j — 1) 
is equivalent to the constraint 

Co 1 = *W ( 5 - 12 ) 



which is interpreted as a cutting off boundary condition for the chain (5.1 ). For simplicity we put 
j = 0, so the boundary condition becomes 

*ri = *o,i- ( 5 - 13 ) 



Remark 1. Hirota equation (5.1) is compatible with the reduction of the type of parity 

tm+i,k = ti,m+k (5-14) 



and boundary condition (5.13) is a consequence of this reduction. 



Following [2T] we can construct a Lax pair for the reduced chain. Under the boundary condition 



(5.13) we have coincidence of the invariants 

Kty(n,m,j) = Ki g (n + l,m,j - 1) 



(5.15) 

K 2l p{n,m,j) = K 2g {n + l,m,j - 1) 

and 

K lg (n,m,j) = K xil (n,m- l,j - 1), 

(5.1b) 

K 2g (n,m,j) = K 2i! (n,m- 1, j - 1) 
and we can relate the eigenfunctions 

<i = A< > 9% = ^lo- (5-17) 



We study the finite reductions of the chain (5.1) on a finite interval < j < Nr. The 



reduction is obtained by imposing the boundary conditions at the left end-point j = Nl 

(5.i8) 
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and respectively at the right end-point j = Nr 



,N R +1 _ ,N R -1 
fc 0,l ~ fc l,0 



(5.19) 



First we concentrate on the left end-point. Due to the reasonings above the eigenfunctions 
should satisfy the following gluing conditions 



(5.20) 



These conditions allow one to close the Lax equations at the left end-point 

N L 

#-1,0 



+N L +1.N L 



+ N L .N L +1 
fc 0,0 t -2,0 



1 L 0,l L -l,l N L 
\ +N L .N L #0,1 • 
t 0,0 t '0,l 



(5.21) 
(5.22) 



From (5.2), (5.6) we have: 

N L 



.N L +1.N L +1 
L 0-l L -l,0 n N L 



#0,-1 — Nl+1 Nl+ i #0,0 
r 0,0 1 



-1,-1 



+ N h +l.N L 

H,0 l 0,0 jN L _ ol N L +l 
+ N L +X.N L 
L 0,0 b l,0 



+N L ,N L +2 

6 0,-l 6 -l,0 N L + l 
+N L + 1.N L + 1 #0,0 > 
t 0,0 -1,-1 



(5.23) 



Vofl ro,0 



To derive similar equations at the point Nr we use the right end-point constraint (5.19), for 
which we have 

,Nr _ „i,N B -l jN R _ N R -1 



i/0,0 ~~ ro,-l > ro,0 ~~ i/l,0 

These conditions allow one to close the Lax equations at the right end-point 



(5.24) 



,N R -l,N R 
J.Nr-I _ l 0,0 r l,0 , N R ~1 _ N R -1 
™1,0 - ,JV«-1,7Vp^0,0 #1,0 > 



H.O ''O.O 



,AT H ,JVr /,iV K -l\2 
..jj-1 _ t -0,-l t '-l,0 Ar fl -1 _ ^0,-1 J /.JVr-I 
#0,-1 — j_N R ,N R #0,0 ,Af n ,AT n ^0,-1 



ATr-1\2 



'UO —1,-1 



J-^R+^R 
L 0,0 



(5.25) 
(5.26) 



From (5.2), (5.6) we have: 



_JV H -1 _ -1,0-1,0 f.JV.R-1 ,JV fl -2\ 

#-i,o - ^atr-i^ath l#o,o "+■ 9o,o ) 



''CO t -2,0 



(5.27) 



26 



<r = <r + 



+N R .N R -2 

L 0,l L 0,0 
j.N R -l.N R - 
L 0,0 L 0,l 



-ih B 
1 ^0,0 



Nr-2 



(5.28) 



The Lax pair found above is not convenient to work with since the operators contain shifts in 
opposite directions. They can be rewritten in the usual form of Zakharov-Shabat system. 



Shift of equation (5.21) forward with respect to the variable n brings it to the form 



+N L +1.N L 
N L _ L -Xfl L 1,Q N L 



01,0 



N L 



N l _iN l+1 00,0 ^0,0 



1 ■ L t 
L 0,0 L 0,0 



(5.29) 



0i,o 



t j+1 t j 

t -l,0 t l,0 j 



+ k 



k=N L J 0,0 r 0,0 

N L + 1 < j <N R -1. 



t 0,0 t '0,0 



Shift of equation (5.26) forward with respect to the variable m brings it to the form 

( + N R -l\2 

l fc o,o 



+ N R ,N R 
N R -1 _ 6 0,1 L -l,0 N R -1 
00,1 — ,N R .N R 00,0 
r 0,0 l -l,l 



,iv R ,i\ R 
L 0,0 L -l,l 



ro,0 > 



(5.30) 



l '-l,0 l '0,l 



00,1 



- j , t '-l,0 t '0,l j+! 

^0,0 '-l,! fc 0.0 t -l,l 



( / -l,0' / 0,l j 

i+T7J+r0o,o+ 2^ 



L a,o L -id 



Nr-1 ,j .k+1 +k+l 
L 0,0 L 



-l,0 t '0,l 



jj'+l + fc +fc+l 
fc=j+l t '-l,l t '0,0 t 0,0 



+N R -l 
k , r 0,0 r 0,0 -JVfl-l 
00,0 ,,-+l,JV n ^0,0 



t 0,0 fc 0,0 



So we have the following system of linear equations 



+1,1 ^0,0 ^0,0 > ly L ^ 7 JV -R z > 



Via 



^O.O t l,0 



(5.31) 



N R -1 



th B 
<Al,0 



j-Nr-I.Nr 

L lfl I Nr-1 



Nr-1 



+Nr-1.Nr ^0,0 
L 1,0 l 0,0 fc=Af L 
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+ E(- 1 )^SA+(-i)^ 



t 0,0 t '0,0 



> 



(5.32) 



, + N l +1.N r -1 Nr-1 ,N L + l.k+l.k+l 

\ L _Y,_ 1 -'0.1. '0.0 , ,-V;v-l^_ 1 111 IU ' -I- >'.,/•• 



^0,1 - ^0,0 + x ,N L ,N R ^0,0 

r o,i r o,o 



^O,! t, 0,0 t U0 



00,0' 



(5.33) 



^^^o + ^ttftCo 1 ' AT L + l< J <iV i? -l, 

*0,0 *0,1 



(5.34) 
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j j-k fk+l 

9lo = £ + (-!) i+1 A<^ ^ < j < N R - 1, (5.35) 

k=N L t O,O t O,0 



j t -l,O t O,l j t -0,O t -l,0 , '0,l fe . t O,O t O,0 /iV H -l , r ^ ■ ^ A r -i /r o«\ 

<i = TJ+ITJ+T^o + 2^ 77+17^1^,0 + Tj+IT^ V>o,o > N L < j < N R — 1. (5.36) 

r 0,0 fc=j+l t -l,l t O,O t O,0 r 0,0 r 0,0 

In the next subsection we gather these equations to a matrix form. 



5.1 Integrable systems corresponding to the Cartan matrices D^, A± 



Imposing of non-degenerate boundary conditions (5.18), (5.19) at both end-points Nl = and 



(2) 

Nr = N leads to the system corresponding to the matrix D N+1 , N > 2 



t, 0,0 ( 'l,l 


t l,0 t 0,l 


= (4,l) 2 , 


f j +3 
i 0,0 i l,l 


t 'l,0''0,l 


- t j ~ l t j+l 


-iN j-N 
( '0,0 t 'l,l 


j-N +N 


= Ko 1 ) 2 



(5.37) 



Write the set of equations (5.31 )-(5.36 ) in the form of a Lax pair for system (5.37). Denote the 

- A _1 ) T . Introduce 2N x 2N matrices 



eigenvector as follow P = (ip°, ip 1 , . . . , ip N 1 , g°, g 1 , . . . , g 1 



( 



A 



fl fV 

t l,0''0,0 
+1 f 
r 0,0 r l,0 





(-l^A 

-A 

A 

(-l) l+1 A 

(-1)"A 







.N +N-1 
r l,0 I 0,0 











N r i,o l: -i,o 
1 -o,o t o,o 



I, 1,0 I '-1,0 



r l,0' : -l,0 
f fl 

1 -o,o t o,o 



-1) 



4.0 fl 

N—l r i,o r -i,o 








j-1 +2 

i \N— 1 t i.o t -i,o 



r 1.0''-l,0 
t 0,CT0,0 



l 0,0''0,0 



I 0,0 C 0,0 

^ -q TV 6 yr- i .. 








\ 



+N-1.N 
r l,0 l -l,0 
.N-l.N 
r 0,0 r 0,0 






fl fZ 



■ iV-l.JV 
1,0 -1,0 
,JV-1,jv 
l 0,0 l 0,0 
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B 



V 



f f 2 

i -o,o'-o,i 

t O,O t O,l 















1 

fl f'i 

t 2 t 2 














1 










1 'o.i'-o.o 

\ fO fN x t° t° t 1 

r 0,l''0,0 t O,0 1 '0,l 0,0 



(A 12 + 1 j-1 j-2 f 

\ l 0,l> 5 -l,0 1 l o,i l -i,o l 








Ti n Wi n 



-1,1 '0,0 



-l,l"0,0 



(Co" 1 ) 2 








0-1,1 







1,0"0,1 

X fO fi fi 

t 0,l , '0,0 , '0,0 








fO f 2 f 2 
'•Q.O -l.C'-P,! 
fl fl f 2 
-1,1 1 '0,0 , '0,0 



-i,o"»,i 



"0,0" 



-1,1 







It is straightforward to check that the compatibility condition of the equations 

P 1)0 = AP, P ,i = BP 



f l j-N j-N \ 

1 r o,i' : -i,o'-o,i \ 

XfO fN-l.N 
t 0,l t 0,0 r 0,0 






t° t 



0,0"-l,0"0,l 



.N-l+N 
l,l l 0,0 b 0,0 



N 



N 



"0,0'--l,0"0,l 
j.2 +N-1.N 
t -l,l''0,0 "-o.o 



jiV j-JV 
""0, 0-1,1 



(5.38) 



is equivalent the system (5.37) 



Example 1 Consider a particular case when Nl = and Nr 



t 0,0 t 'l,l 



t 1 t 1 



t, l,O fc 0,l 

t 1 



T/ie system corresponds to algebra A\ 



(i) 



(5.39) 



1,0"0,0 



.4 



0,0"l,0 

-A 



+ A 



r l, 0-1,0 



fU fl 



*0 fl 

t 0,0 t 0,0 



/fo Lax pazr is given by (5. 38) with the matrices 
1 

B 



1 1 r (),l 0,0 


1 


(*o,i) 2 *-i,o 


A t 0,l t 0,0 


A 


+0 j-U fl 
I 0,0 I 0,1 !: 0,0 






'0,1 -1,0 


""0, 0-1,1 




l 0, 0-1,1 



Example 2 // one imposes condition (5. IS) for Nl 



t 0,0 , 'l,l t l,0 , '0,l 



t 'Q,O fc l,l t, l,O t '0,l 
t 0,0 fc l,l L 1,0 L QA 



and ( f5.igp /or N R = 2, then 
(4,i) 2 , 

t l,0''0,l> 



(5.40) 



The Lax pair is of the form (5.38) where A and B are 4x4 matrices 
( 



A 



fl t o 

t 0,0 t l,0 



-A 
-A 



-1 



"1,0"0,0 
f2 fl 
I 0,0 I 1,0 







"1,0"-1,0 
r 0,0 r 0,0 







\ 



t 



0,0 ! '0,0 



A 








"1,0-1,0 
*0 +1 
0,0 0,0 
fO fl 
n, 0-1,0 

t 0,0 t 0,0 



r 0,0 r ' 





ft 



fl 



B 



fl fl , f+1 ^2+1 . +1 f2 f 2 \ 

1 '-0,l'-(),0 1 V-0,1- 1 -1,0 1 '■Q.l'-l.o'-O,! ' 



AtS,i 



+i ti 

I 0,0 r 0,l 



0,1"0,0 

1 








t fl 



X fO fO fi 

A t 0,0 t 0,l l '0,0 



A t 1 i 2 
t 0,l 1, 0,0 , '0,0 



"0,0"0,0 



t 2 

t -l,l''0,0 

(^,o) 2 

-l.l'O.O 



t 0,0 t -l,l 





0,O t, -l,0'-0,l 



i 1 



il 

i.i'o.o'-o.o 



f 2 f 2 
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5.2 Integrable systems corresponding to the Cartan matrices A 



N 



In addition to boundary condition (5.13) we introduce the degenerate boundary conditions 
at the left end-point and 

t N R +l = q t N R = j 



(5.41) 



(5.42) 

at the right end-point. The degenerate boundary conditions imply that the corresponding eigen- 
functions are zero: g NL ~ Y = and i/j Nr = 0. 

In order to obtain the system corresponding to the Cartan matrix An-i we cut off the Hirota 



chain by imposing degenerate boundary conditions (5.41) at the point Ni = and (5.42) at 
Nr = N. The resulting reduction is as follows 



t 1 t 1 

t, l,O t Ul 



t 0,0''l,l ''l.C^Ojl 

j.JV-Ij.JV-1 _ j.JV-Ij.JV-1 
l 0,0 c l,l c l,0 C 0,l 



t 2 



t{; 1 C 1 1 ,2<j<N-2, 



(5.43) 



WV-2 

r i,o • 



In this case our algorithm gives the Lax pair found years ago in [13]. In order to formulate it 
introduce the eigenvector P = (-0°, ip 1 , . . . , if; N ~ l ) T and N x N matrices 

\ / 



*0,0 





-1 

.2 .1 

r l,0 I '0,0 
72—7]— 
r 0,0 r l,0 





V 






















-1 










t l,0 I 0,0 
t 0,0 l l,0 


-1 















t N-l 

.N-l 
c l,0 



V 







l 0,l 



1 t 1 

0,0 t 0,l 



t 1 t 




t 0,l t, 0,0 

t 2 t 2 






1 



V 

















,N-2 
I 0,0 



j iv — i j i\ — 
c 0,0 r 0,l 






1 

(5.44) 



It is easy to check that the compatibility condition of the equations 

P 1)0 = UP, P ,x = VP 



(5.45) 



leads to the system (5.43) 



30 



5.3 Integrable systems corresponding to the Cartan matrices 



We impose boundary condition (5.41 ) for Nl = and (5.19) for Nr = N. The resulting reduction 
is as follows 



t 1 t 1 

{, 0,O t l,l 


-t 1 t 1 

t l,0 t 0,l 


- t 2 




fj fj 

t 0,0 t 'l,l 


-fj fj 


= 1&4f,2<j<N-l, 


(5.46) 


j.N j-N 
t 0,0 t 'l,l 


j-N 4-N 
£ 'l,0 t 0,l 


= (fjo" 1 ) 2 . 





System (5.46) can be rewritten in form of (1.3) by changing the variables P ' — e u . It corresponds 
to the algebra B N . The system admits a Lax pair. Let us denote P = (ip°, • • • > 4' N ~ 1 : 9°> 9 X ■, ■ 



g N- T 



and introduce 2N x 2N matrices 

*0,0 





-1 

.2 .1 
t l,0 t 0,0 













f t 1 

l '0,0''l,0 












-fiV j.N-1 
t l,0 t 0,0 



o.on.o 








-1) 

t 1 



N z -i,o 
fl 
l 0,0 



■1) 






il +2 

N—l t i,o i; -i,o 
t 1 t 2 

t 0,0 t 0,0 



fl 



fl 

l 0,0 



'-l.cr-i.o 

t 1 t 2 
t o,o t o,o 



V 

/ 



' i I 



t 0,0 t 0,l 





V 





o 
o 

o 







1 

+1 t 3 
t 0,0 1 '0,l 
72—72— 
r 0,0 r 0,l 












1 



,N-2,N 
"-0,0 "-Q,! 
.N-l.N-l 
0,0 r 0,l 













'•-i.i'-o.o 

.1 ,N-1 

t 0,0''0,0 

I 2 t N 

t -l,l''0,0 



v r o,o > 



0.1-1,0 
■•0, 0-1,1 





1 (1 J.Z 

-l,l r 0,0 r 0,0 

t 2 i 2 

t -l,0 t 0,l 
















+N-1+N 
hjO -1,0 
j-N-l.N 
l 0,0 l 0,0 









-1*-1,0 

+1 

l 0,0 


(- 


1)* 


f l f 2 

H., 0-1,0 
f l f 2 
t 0,0 I '0,0 


.N-l.N 

n,o t -i,o 

t N-l,N 
l 0,0 l 0,0 












































1 














j-1 j-N-l.N 
t -l,l''0,0 l 0,0 
j-1 j-N j-N 
L 0, 0-1, 0^0, 1 
j.2 j-N-l.N 
-l.l'-O.O 0,0 



Then according to our general scheme the compatibility condition of the equations 

P 1)Q = AP, P 0)1 = BP. 



(5.47) 



(5.45 



(5.49) 
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leads to the system (5.46). 

Remark 2. The system B>n can be obtained from the system A 2 n-i by imposing the cutting 
off constraint of the form 

= (5.50) 

(see Remark 1 above). 

6 Method of finding integrals from Lax representation for 
the systems corresponding to the Cartan matrices An, 

B N 



For any integer N > 2 system (5.43), corresponding An, admits a complete set of independent 
integrals consisting of N independent m— integrals and N independent n— integrals. Demonstrate 
an algorithm of finding them by using Zakharov-Shabat operators. Begin with m— integrals. Due 
to the definition we have two different expression for the shifted eigenvector Pk,x, k > 1 

Pk,l = VkflUk-lflUk-2,0 ■ ■ ■ UiflUofiP, 

and similarly 

Pk,i = £4-i,if4-2,i • • • t4,i£4,iVo,o-P- 
Comparison of these two formulas yields 

-Dm (^4- 1,0 ^4- 2,0 • • • £4,0^0,o) = Vfc,0^4-l,oC4-2,0 • • • ^1,0^0,0 ^o.o 1 - (6- 1 ) 

Due to the triangular structure of the matrices V^o and Vo,o the map converting any upper trian- 
gular matrix X to a matrix X = V^^XV^q leaves unchanged the element of the matrix X located 
at the right upper corner: (X)i t N = -^i,jv- Thus the corresponding element, denote it through 
Ik, of the upper triangular matrix £4-1,0^4-2,0 • • ■ ^1,0^0,0 is an m-integral. It can be proved that 
integrals In, In+i, I2N constitute a complete set of m— integrals. In a similar way one can find 
integrals in the other direction. 

Let us illustrate the statement above with the following 
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Example 3 Consider the system (5.45) for N = 2 



Recall its Lax pair 



where P = (if) , ip ,ip ) and 



( * 



U 



\ 



L 0,0 








t, 0,0 I 'l,l ,/ l,o''0,l 



f 2 f 2 
t 0,0 t 'l,l 



t l,O fc O,l 



"0,1' 

+1 



UP, P ,i = VP, 



t l,0 ! '0,0 
c 0,0 r l,0 





"-0,0 



V 



' n 



V 



t 0,0 I '0,l 





l 0,0 



t O,O t O,l 



0^ 




(6.2) 



(6.3) 



(6.4) 



Evaluate the elements at the right upper corner for the following two products ^2,0^1,0^0,0 ond 
^3,0^2,0^1,0^0,0 an d find two independent m-integrals 



P 



^3 _|_ ^2 _|_ ^0 
t 2 ^2^1 ^1 



l + l/ + 2\2 



"4 7 1 ^ 
1 W+l 



*ai*i(* 



■ , 1 



(/ 3 (/ 2' / _|_ ^0 

(t 2 ) 2 t 3 t 1 ^3^1 ^2 



where the second index for the variables t ,t is omitted 



In a similar way integrals for the system (5.46), corresponding to the algebra B^, are constructed. 



Let us consider a matrix $ of dimension 2N x 2iV: 




(6.5) 



where En is the unity matrix of dimension N, P12, -E21 are matrices of dimension iV with zero 
entries, P22 is a matrix of dimension iV with the unity adverse diagonal. By an automorphism 
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X — > $X$ the matrices (5.47) and (5.48) are transformed to the triangular matrices 



A 



( |o 

*0,0 











B 



"0,0"0,1 



V 





-1 





o 





t i,o' , o,o 

,N ,N-1 

r 0,0 r l,0 






,JV-1,jv 

n.o - 













,iV-l,JV 
'•o.o c 0,0 
+N-1+N 

n.o ■'—l.o 



.N-2.N-1 

HjO z -i,o 

.JV-2.JV-1 
r 0,0 I 0,0 

.N-2.N- 

fi n £ 



1 

1.(1 



iV-1 .JV 
0,0 c 0,0 



,JV-2 ,]V-1 



-1,0 





1 








.. 

.. 

.. 

.. 



,N-2,N 

.N-l.N-l 
l 0,0 l 0,l 













1 

.N-l 
l 0,0 













-l.l'O.O 



JV-1n2 



(*0,0 



,N-2,N j-N 
z 0,0 r -l,0 r 0.1 
AT-l .JV-1.JV 



.N-l+N 
l.l'-O.O l 0,0 



,JV-2 ,JV-1 

-i.i'o.o "-o.o 








*0,0 
AT-l f -i,0 

+1 

l 0,0 

■1)^4^ 

*0,0 



\ 



-1) 



<-o,o 












t 1 t 1 

l 0, 0-1,1 



/ 
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(6.6) 



By an automorphism acting as follows X — > F 1 XF matrix (6.6) is transformed to a lower 



triangular matrix for which all diagonal entries are equal to the unity. Here F is a 2N x 2N matrix 
of the form 

En E12 
E21 F22 



(6.7) 



F22 is a diagonal matrix such that F22 = diag . - 

V -1,0 * 



/ ,N ,N-1 ,1 
/ r 0,0 r 0,0 I 0,0 
I t N ; ,N-1 ) • • • , II 



and are defined in (6.5). 



Reasonings similar to that of the case A^r allow one to derive the integrals. 



Example 4 Consider the system (5.46) for N = 2 



/i +1 _ +1 +1 _ + 



t 0,0''l,l ''l.CrO,! 



0,1' 



(t{ Y 



(6. 
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Recall its Lax pair 



P 1>0 = AP, P , 



BP, 



(6.9) 



where P = , ip l , g° , g 



1\T 



and 



A 



( *i.o 
71 
c o.o 







.2 .1 

r l,0 I 0,0 



t O,O l l,0 






+1 
'().() 

fl 

-i,o 
71 

l 0,0 



71 

r 0.0 



fl f2 



fl f 2 



7!— 72 

c O,0 I 0,0 



B 



l O,(rO,l 










1 



fl 

l O,0 








t 1 t 1 

■•0,1 - 



1,0 



-l,l'-0,0 

(*o,o) 2 



71 +1 

c 0, 0-1,1 








t 2 

-i,(ro,i 



1 fl y.2 

-1,1 o,o' : o,o 

.2 ,2 
-1,0 C 0,1 
72—72 



(6.10) 



an automorphism X — >■ $X$ transform the matrices A and B to the triangular matrices 



A 



( io 

+1 
'o.o 







1 







f 1 t 

i n L 



1,0'-1,0 



f2 fl 



f 1 t 2 
t o,o'-o,o 



V 








t 



0,0 r 0,0 









+1 

l 0,0 



"0,0 



/ 



B 



J 



tl 

0,1 



f 1 f 1 





V 









1 

(^,o) 2 

72 72— 
--1,1 C 0,0 

+i 
r o,o 

t l f2 








,2 ,2 
-l,0 r 0,l 



-l,l t 0,0 t 0,0 









t 1 t 1 
-i.i'o.o 



/ 



(6.11) 



By automorphism X — > F l XF transform the matrices (6.11) to the triangular matrices 

( *i.o 



.4 



V 



71 

r o,o 







fA Jl 







t 0,0 l l,0 






fl 

f l f 2 





1 



\ 



fl 



f2 

1 t'l 



I- o,o'"i,o 





t 2 

l l,0 



/ 



B 



J 



1 



t 1 t 1 

t 0,0 t 0,l 



V 










1 

(*o,o) 2 



f 2 f 2 
l 0, 0^0,1 

fl 

r 0,0 






1 



fl 



l o,i'-o,o 



f 1 f 1 



(6.12) 



/ 



Evaluate the elements at the right upper corner for the following products ^2,0^1,0^0,0 an d 
-^3,0-^2,0-^1,0-^0,0 an d find two independent m-integrals 



I 1 = -4- 



tl 



fl/2 



f l f 2 



t 1 ' 



tl 


t\t 2 t^ 


f l f 2 
L l u 3 


f l f 2 f l 

L 2 L L 4: 


(tl^tytl 




tl 


t\t\t\ 


fl/2 
L 3 L 1 


t\t\t\ 


t\t\t\t2 






t l t 2 


tH 1 


t l t l t 2 


f l f l f 2 


th 




f l f 2 
L 1 L 2 


tH 1 


t 1 t 1 t 3 


t 1 t 1 t 2 


tl 



Since the integral I 2 is too complicated one can choose a more simple one 



P = i 2 + Pdj 1 = § + 



t 2 



(^i) 2 ^2 _|_ 2_^ML _|_ 



(mi t 



2^+2 
l 2 



(6.13) 



(6.14) 



(mi (tiy (mi 

In the integrals above the second index for the variables t x ,t 2 is omitted, because its values is zero 
for all considered variables. 
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7 Periodic boundary conditions 



In this section we discuss briefly the well known periodically closed reduction of the Hirota chain 
(see for more details and the references 
boundary conditions 



f 



t 



N 



. Close the chain (5.1) by imposing the periodical 

(7.1) 



t 



N+l 



t°. 



Close the Lax pair (5.2) by setting the conditions on the eigenfunction 

1 



i)- x = Xi/j 



N 



■0 



As a result we get a finite system of the form 



+0 +0 _ +0 j-0 _ +N +1 

t '0,0 t 'l,l b i,o L o,i — tl nt 



N+l 



A 



(7.2) 



1,0^0,1' 



t 0,0 , 'l,l t 'l,0''0,l ~~ ^1,0 ^O,! > 1 — J — ly x ' 



(7.3) 



+ Af _ + N ±N _ + j-1 
'UOn,! t l,0 t '0,l ~~ t l,0 , '0,l 



which is closely connected with the Cartan matrix A N . 

Boundary conditions (7.2) reduce the sequence of linear discrete equations ( |5.2 ) to a Lax pair 
for the reduced system (7.3 ). Introduce the eigenvector P = xjj 1 , . . . , ip N ) T and N + 1 x iV + 1 
matrices 

\ / 



U 



/ t° t 1 

/ r O.0 r l,0 
.0 .1 

r l,0 c 0,0 






-1 



r l,0 r 0,0 
t 0,0 I, l,0 







-1 

f 3 f 2 



V 



o 











'•o.cri.o 

f N ,0 I 



V 



1 



r 0.1 r 0,0 
t 0,0 , '0,l 







/d +1 





1 



V 



o 















f N-l f o 



j-N f N 
t 0,0 r 0,l 



.1 ,N \ 
\ r 0,l 0,0 \ 

c o,o t o,i 




1 / 



If the functions t°, t 1 , . . . , t solve the system (7.3) then the overdetermined system of linear 
equations 

P 1>0 = UP, P ,i = VP (7.4) 

is compatible (see [32]). 



Remark 3. Note that system (7.3) slightly differs from the system corresponding to the same 



algebra Aw but defined by the formula (1 1 . 3h . For example the system (17.31) for iV = 2 is of the 
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form 



t 0,0 t 'l,l 



In terms of the variables u c 



t 1 t 1 

£ '0,O t 'l,l 

t 2 t 2 
logt , u l 



t 1 t 1 



t 2 t 2 



t 2 t 1 
L l,0 L 0,n 

f\ 4.0 

t l,0 t 0,l' 



(7.5) 



— logt 1 , u 2 



logt it looks like 



A(u°) 
A(u x ) 
A(u 2 ) 



g u l,0+ u O,l~ 



<4,r 



4 i,o 



(7.6) 



1,0" 



'•*0,1 



while formula (1.3) gives the system 



A(u°) 
A(u x ) 
A(u 2 ) 



4, - 



l l,0~ 



'4,i- 



4,1 



(7.7) 



u l,0 U 1,0+ U 1, 0+^0,1 



Let us give also the Lax pair for the system (7.5). Let us consider vector P = , ip 1 , ijj 2 ) T . 
Introduce 3x3 matrices 



/ t° t 1 

f fl 
t l,0 t 0,0 



u 



-1 

+2 fl 
t l,0 t 0,0 
72—71— 
r 0,0 r l,0 







\ 



V 



c 0,0 r l,0 



/ 



V 



J 



1 



t 0,l t 0,0 



fl fl 



V 





1 

fO f l 
r 0,l I 0,0 

73—73— 

I 0,0 t 0,l 



A 



fl fZ 

t 0,0 t 0,l 



1 



\ 



If the functions t°,t l ,t 2 satisfy the system (|7.5|) then the equations 



Pi,o = UP, P ,i = VP. 



{71 



are compatible. 



Remark 4. The Lax pair (7.4) can be rewritten in terms of the Cartan-Weyl basis of the 



algebra Al 1 ^ 



N ■ 



P 10 = (-A + e tfi,o-Ub,o) p f 

p 0>1 = (E + e Uo >iAe- Uo '°) P 



(7.9) 



N _ N N 

where A = £ fa A = £ U = £ e""* = f 

i=0 i=0 i=0 



= 0, [ei,/j] = Sijhi 



(7.10) 
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[hi, ej] — AijCj, [hi, fj] — —Aijfj, (7-H) 



and Aij are elements of the Cartan matrix of the algebra A$. 



8 Conclusions 

A map is suggested converting any N x N matrix A to a finite system of difference-difference 



equations of exponential type (see (1.3)). A hypothesis is formulated claiming that if A coincides 
with the Cartan matrix of a finite or afline Lie algebra then the corresponding system is integrable. 
The hypothesis is approved by numerous examples. The systems obtained are rather simple and 
elegant. They essentially differ from those studied earlier (see survey [XSJ and references therein). 
For instance, the system corresponding the algebra G2 given in [16] reads as: 

t£\u-1)t£\u + 1) = T^Uu)T^Uu)+T^(u), 
T?2 (u - I) ?2 (« + I) = T^.MT^u) + T« (u - 1) T£\u)T^ («+§), 
T ?2 + i (« - \) T 3 2 2 + i (« + §) = ^(«)T& a («)+r,k 1) (tt-i)TW(ri+i)ia. 1 (tt), 

^3m+2 { U ~ |) ^3m+2 ( U + |) = ^3m+l ( M )^3m+3 ( U ) + ^«^( n )^~m+l ( M ~~ |) ( n + |) ) 

(8.1) 



while our formula (1.3) generates G2 system which can be represented as follows: 



t 1 t 1 -t 1 t 1 = t 

v n,m l/ n-\-l,m-\-l n+l,m"n,m+l n,m+l' 

t 2 t 2 -t 2 t 2 = (t 1 ) 3 
n,m n+l,m+l n+l,m n,m+l V n+l.m/ 



(8.2) 

n,m n+l,m+l n+l,m n,m+l 



In a recent article by Kimura, Yamashita and Nakamura [31] a new application of conserved 
quantities of discrete-time integrable systems to numerical computations is suggested. The systems 
studied in the present paper might have applications in such kind numerical methods. 
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